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Coordinate transformations can be used to manipulate fields in a variety 
of ways for the Maxwell and Helmholtz equations. In Sect. [T] we focus on 
transformation elastodynamics. The idea is to manipulate waves in an elastic 
medium by designing appropriate transformations of the coordinates and 
the displacements. As opposed to the Maxwell and Helmholtz equations, the 
elastodynamic equations are not invariant under these transformations. Here 
we recall the transformed elastodynamic equations, and then move to the 
effect of space transformations on a mass-spring network model. In order 
to realize the transformed networks we introduce "torque springs" , which 
are springs with a force proportional to the displacement in a direction other 
than the direction dictated by the spring terminals. We discuss some possible 
homogenizations of transformed networks that could have applications to 
manipulating waves in an elastic medium for e.g. cloaking. 

Then we look at an approach to cloaking which is based on cancelling 
the incident field using active devices (rather than passive composite materi- 
als) which are exterior to the cloaked region. Exterior means that the cloaked 
region is not completely surrounded by the cloak, as is the case in most trans- 
formation based methods. We present here active exterior cloaking methods 
for both the Laplace equation in dimension two (Sect. [2]) and the Helmholtz 
equation in dimension three (Sect.[3|. 

The cloaking method for the Laplace equation we present in Sect. [2] applies 
also to the quasi-static (low frequency) regime and was in part presented in 
[251 We first reformulate the problem of designing an active cloaking 
device as the classic problem of approximating analytic functions with poly- 
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nomials. This theoretical approach shows that it is possible to cloak an object 
from an incident field with one single exterior device. Then we give an ex- 
plicit solution to the problem in terms of a polynomial and determine its 
convergence region as the degree of the polynomial increases. This conver- 
gence region limits the size of the cloaked region, and for the new solution we 
propose here it allows one to cloak larger objects at a fixed distance from the 
device compared to the explicit polynomial solution given in [231 IE] • We also 
discuss how our approach can be modified to simultaneously hide an object 
and give the illusion of another object, in the same spirit as illusion optics 

Next in Sect. 3] we consider the Helmholtz equation and use the same 
techniques as in [22j to show that in dimension three it is possible to cloak 
an object using four devices and yet leaving the object connected with the 
exterior. Our method is based on Green's formula, which ensures that an 
analytic field can be reproduced inside a volume by a carefully chosen single 
and double layer potential at the surface of the volume. Then we use addition 
theorems for spherical outgoing waves to concentrate the single and double 
layer potential at a few multipolar sources (cloaking devices) located outside 
the cloaked region. We determine the convergence region of the device's field 
and include an explicit geometric construction of a cloak with four devices. 

The three sections of this chapter can be read essentially independently of 
each other. 



1 Transformation elastodynamics 

Transformation based cloaking was first discovered by Greenleaf, Lassas and 
Uhlmann |151I16) in the context of the conductivity equations. Independently, 
Leonhardt realized that transformation based cloaking applies to geometric 
optics [2H] and Pendry, Schurig and Smith [l^ realized that transformation 
based cloaking applies to Maxwell's equations at fixed frequency, and this 
led to an explosion of interest in the field. It was found that transforma- 
tion based cloaking also applies to acoustics [HI El IIZ] , which is governed by 
the Helmholtz equation, provided one permits anisotropic density |45| . These 
developments, reviewed in [H [TSl HI jB] rely on the invariance of the con- 
ductivity equations. Maxwell's equations, and the Helmholtz equation under 
coordinate transformations, and have been substantiated by rigorous proofs 
[171 EH ES] . The invariance of Maxwell's equations under coordinate trans- 
formations has led to other envisaged applications such as field concentrators 
[H], field rotators [7], lenses [15], superscatterers [SD] (see also [33]) and the 
name "transformation optics" is now used to describe this research: see, for 
example, the special issue in the New Journal of Physics [30] devoted to 
cloaking and transformation optics. The perfect lens of Pendry [42] can be 
viewed as the result of using a transformation which unfolds space [29) and 
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associated with such folding transformations is cloaking due to anomalous 
resonance [5^ 1 HO I [55). 

A largely open question is how to construct metamaterials with the re- 
quired combination of anisotropic electrical permittivity £(x) and anisotopic 
magnetic permeability /^(x) needed in transformation optics designs, fre- 
quently with £(x) = /x(x). Only recently was it shown [34], building upon 
work of Bouchitte and Schweizer [2], that any combination of real tensors 
(e, fi) is approximately realizable, at least in theory. 

Curiously, the usual elastodynamic equations do not generally keep their 
form under coordinate transformations. Either new terms enter the equations 
[37] , so they take the form of equations Willis introduced [49] to describe the 
ensemble averaged elastodynamic behavior of composite materials (which are 
the analog of the bianisotropic equations of electromagnetism [47j ) . or the 
elasticity tensor field does not retain its minor symmetries [3|. Nevertheless, 
as shown in and as is explored further here, there is some hope that meta- 
materials can be constructed with a response corresponding approximately 
with that required by the new equations. 



1.1 Continuous transformation elastodynamics 

By extending the analysis of [37], let us show that the equation of elastody- 
namics 

- V • (C(x)Vu) = a;2p(x)u (1) 
changes under the transformation 

x' = x'(x), u'(x'(x)) = (B^(x))-iu(x) (2) 

to the equation 

- V • (C'(x')V'u' + S'(x')u') + D'(x')V'u' - cj2(p'(x')u') = (3) 

where the tensors C, S', D', p' are given in terms of the functions x', B and 
their derivatives. Here the transformation of the displacement is governed 
by B(x) which can be chosen to be any invertible matrix valued function. 
(The inverse and transpose in (B"'"(x))~^ have been introduced to simplify 
subsequent formulae.) 

Indeed let us first note that 

dxi dxi dxi dx'^ dxi ^ 

= A'^(V'u')B + G'u' (4) 

in which A and G are the tensors with elements 
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dx' dB 
" dx, ' ~ dx, ■ ^ ' 

Now ([T]) implies that for all smooth vector- valued test functions v(x) with 
compact support in a domain i7, 

= / hV • (C(x)Vu) - wV(x)u] • V dx 

= / [C(x)Vu : Vv - w^p(x)u • v] dx 

= / [C(x)(A^(V'u')B + Gu') : (A^(VV)B + Gv') - a;2^(x)(B^u') • (B^v')]a"^ dx' 
= / [C'(x')V'u' : VV + S'(x')u' : V'v' + (D'(x')V'u') • v' - Lo^p' {x.')u') ■ v'] dx' 
=^ [ [-V • (C'(x')V'u' + S'(x')u') + D'(x')V'u' - a;2(p'(x')u')] • v' dx' 

JO' 



(6) 



in which the test function v(x) has been transformed, similarly to u(x), to 

v'(x'(x)) = (B^(x))-V(x), (7) 

and a(x'(x)) = detA(x) while C'(x'), S'(x'), D'(x') and p'(x') are the 
tensors with elements 



^kij ~ ^ GpqkAirBjgCpqrs — ^ijkJ 

p'ij = O, ^Bii^Bjkp—a ^UJ '^GpqiGrsjCpqrs 



a ^B,kBjkp-a ~Q^~Q^^vqrs- (8) 



,dB^dB 

From (|6]) we see directly that ([T]) transforms to (|3|. 

Remark 1. The transformed elastodynamic equation ([3| can be written in the 
equivalent form of Willis-type equations [JS] 

V' • cr' = -iwp', 

cr' = C'(x')V'u' + (^/w)S'(x')(-^LJu'), 

p' = p'(x')(-W) + (^/w)D'(x')Vu', (9) 

in which the stress cr', which is not necessarily symmetric, depends not only 
upon the dispacement gradient V'u' but also upon the velocity —iujxi', and 
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the momentum p' depends not only upon the velocity — iwu', but also on the 
displacement gradient V'u'. 

Remark 2. If we desire the transformed elasticity tensor C'(x') to have all 
the usual symmetries of elasticity tensors, namely that 

^"ijke = ^'jike = ^'k£ij^ (10) 

then we need to restrict the transformations to those with B = A. This was 
the case analyzed by Milton, Briane and Willis [37] , 

Remark 3. In the particular case where B = I the transformation (|8| reduces 
to 

C[^^, = a-^A,pAkrCp,ri, S' = D' = 0, p' = a'^pl, (11) 

corresponding to normal elastodynamics, with an isotropic density matrix p' , 
but with an elasticity tensor C only satisfying the major symmetry C^' ^.^ = 
^'kiij- This was the case analysed by Brun, Guenneau and Movchan [3] in a 
particular two-dimensional example. 

Having derived the rules of transformation elasticity, one can then apply 
the same variety of transformations as used in transformation optics, includ- 
ing cloaking and folding transformations. The point is that a wave propa- 
gating classically in the classical medium can have a strange behavior in the 
new abstract coordinate system x'. If we are able to design a real medium 
following a system of equations equivalent to the transformed system, then 
we are able to force a strange behavior for waves in real physical space. 



1.2 Discrete transformation elastodynamics 



There is a discrete version of the transformation (11). Suppose we have 
a network of springs, possibly a lattice infinite in extent, with a count- 
able number of nodes at positions xi, X2, X3, . . . , x„ . . ., at which there 
are masses Mi , M2 , My, , . . . , M^ . . . , and at which the displacements are 
Ui, U2, U3, . . . , u„ . . .. Let kij denote the spring constant of the spring con- 
necting node i to node j. There is no loss of generality in assuming that all 
pairs of nodes are joined by a spring, taking kij = if there is no real spring 
joining node i and j. Let ^ denote the force which the spring joining nodes 
i and j exerts on node i. Hooke's law implies 

Fij = -Fjvt = fcijnij[nij • (u^ - u^)], (12) 



where 
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is the unit vector in the direction of Xj — Xj. In the absence of any forces 
acting on the nodes, apart from inertial forces, Newton's second law imphes 



-AWu,. (14) 



Now let us consider a transformation x' = x'(x) with an associated inverse 
transformation x — x(x'). Under this transformation the position of the 
nodes transform to x'j^, Xj, Xg, . . . , x^ . . ., where x^ — x'(xi). We focus, for 
simplicity, on the case corresponding to B = I where the forces, masses and 
displacements transform according to 

F',^^. =F,,j, M; = Af„ u^ = u,. (15) 

After the transformation, Newton's second law clearly keeps its form, 

^F^,^. = -M;c.2u^, (16) 



while ( 12 1 transforms to 

F;,^ = -F;,, = fc:,^v:^^K^.(u;-u^)] (17) 



where 



, _ x(x;.)-x(xO 

K3 - - |x(x')-x(xO| 



Hence in the new coordinates x^ the system is governed by equations 
similar to the classical system of equations for a network of masses joined 
by springs, but the response of the springs does not anymore correspond 
to normal springs. While the action-reaction principle F^- ^ — — F^ ^ remains 
valid, the force F^ ^ is not generally parallel to the line joining x^- with x,^. 

Now we desire to construct a real network having a behavior governed 



at a fixed frequency, by the system of equations (16 1 and (17). To that aim 



we need to construct a two-terminal network made of classical masses and 



springs which has the response (17) for any unit vector v.^^ . We call these 
two-terminal networks "torque springs" since they extert a torque in addition 
to the usual spring force. We show how they can be constructed for fixed 
frequency w in the next section. 



1.3 Torque springs 



A torque spring, being a two-terminal network with a response of the type 
(17), is characterized by two terminal nodes Xi, X2, the direction of exerted 
forces vi 2 which can be different from the direction of the line joining Xi 
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and X2 and the constant of the spring ki^2- The existence of torque springs 
is guaranteed by the work of Mihon and Seppecher [35] which provides a 
complete characterization of the response of multiterminal mass-spring net- 
works at a single frequency. The complete characterization of the response 
of multiterminal mass-spring networks as a function of frequency was sub- 
sequently obtained by Guevara Vasquez, Milton, and Onofrei [21]. Here we 
are just interested in constructing two terminal networks with the response 
of a torque spring. In this case a simpler construction, than provided by the 
previous work, is possible. 

Consider the network of Fig. [T] For its design we start with xi, X2 and a 
unit vector v = vi2 not parallel to xi — X2. (A normal spring can be used if v 
is parallel to Xi — X2.) Choose p > and define yi = xi + pv, y2 = X2 + pv, 
and choose a vector w 7^ in a direction different from v and X2 — xi . Define 
zi = yi + w, Z2 = y2 + w, ti = zi + V, t2 = Z2 + V. The pairs (xi,yi), 
(X2,y2), (yi,y2), (yi,zi), (y2,Z2), (zi,Z2), (zi,ti), (z2,t2) are joined with 
normal springs of constant k. Masses (with mass m, where the lower case m 
is used to identify them as internal masses of torque springs) are attached to 
the nodes ti and t2 only. All nodes but Xi and X2 are interior nodes which 
means that no external forces are exerted on them. 

Let us denote by T the tension in the spring (xi, yi), taken to be positive 
if the spring is under extension and negative if it is under compression, i.e. 
the spring exerts a force vT on the terminal at xi and a force — vT on the 
node at yi. Then the balance of forces at node yi fixes the tensions T' , T" in 
the springs (yi,y2), (yi,Zi) in a purely geometrical way. It is easy to check 
that the balance of forces at y2, Zi, and Z2 gives tensions —T, — T", —T', T, 
-T in the springs (x2,y2), (y2,Z2), (zi,Z2), (zi,ti), (z2,t2) respectively. 

All the tensions being determined when one is known, the truss is rank 
one : there is only one scalar linear combination of the displacements Ui, 
U2, Wi, W2 of nodes Xi, X2, ti, ±2 which influences T, and T = if and 
only if this scalar linear combination vanishes. It is easy to check that this 
combination is (u2 — Ui — W2 -|- Wi) • v since displacements leaving this zero 
(floppy modes) do not produce any tension in the springs, as they leave the 
spring lengths invariant to first order in the displacements. Hence there exists 
a constant K (proportional to k) such that T = K{u2 — Ui — v^2 + wi) • v. 
Finally Newton's law ( |l4| gives at nodes ti, t2 respectively T = mw^Wi • v 
and —T = ma;^W2 • v and so T = K{u2 — Ui) • v -I- 2TKm~^uj~^ from which 
we conclude that 

_ KrauP' , ^ , , 
^^U2-Ui -v. 19 

The forces Fi and F2 which this torque spring exerts on terminals 1 and 2, 
respectively, are therefore 

Fi = -F2=rv = fc'v[v(u2-Ui)], with A:' = — , (20) 
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which is exactly of the required form ( |17[ ). If we want k' to be positive then 
we should choose m and K so that muT — 2K > 0. 

There are many other constructions which produce torque springs. An- 
other configuration, which is closer in design to a normal spring, is that given 
in Fig. [2] In two-dimensions this type of construction may be preferable to 
that in Fig. [l] to reduce the number of spring intersections when assembing 
a network of torque springs. It also may be preferable if we wish to attach a 
torque spring say between two parallel interfaces. 

The torque springs described here are quite floppy. To give them some 
structural integrity one would need to add a scaffolding of additional springs, 
extending out of the plane if the torque springs are going to be used in a 
three dimensional network. Provided the spring constants of these additional 
springs are sufficiently small, this can be done with only a small perturbation 
to the response of the torque spring, as shown in [3T] . 

In assembling a network of torque springs it may happen that an interior 
spring or interior node of one torque spring intersects with an interior spring 
or interior node or interior node of another torque spring. Since we have the 
flexibility to move the interior nodes of each torque spring we only need be 
concerned with the intersection of two springs, or between the intersection of 
one spring and a node. In three dimensions if a spring intersects with another 
spring or a node we can replace one or both springs by an equivalent truss of 
springs to avoid this situation. In two dimensions if a spring intersects with 
a node we can again replace the spring by an equivalent truss to avoid this 
situation. Then if two springs intersect in two dimensions they must either 
overlap or cross: if they overlap we can replace each by an equivalent truss of 
springs, while if they cross we can (within the framework of linear elasticity) 
place a node at the intersection point and appropriately choose the spring 
constants of the joining springs so that they respond like two non-interacting 
springs - see example 3.15 in Milton and Seppecher |36) . 




Fig. 1 Sketch of a torque spring. The open circles represent terminal nodes, and the closed 
circles could be either terminal nodes or interior nodes with masses attached. The straight 
lines represent springs. The large arrows represent external or inertial forces acting on the 
nodes at one instant in time. The two small arrows on each spring give the direction of the 
force which the spring exerts on the node nearest to the arrow. 
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Fig. 2 An alternative construction of a torque spring. The straight lines represent springs 
and the circles, large arrows, and small arrows have the same meaning as in figure 1. 



1.4 Homogenization of a discrete network of 
torque springs 



As shown in Sect. |1.2[ the original network of springs with nodes at posi- 
tions Xi, X2, X3, . . . , x„ . . . and spring constants kij responds in an equiva- 
lent manner to the new network of torque springs with nodes at positions 



x']^ , X2 , X3 , . . . , x^ . . . and torque spring parameters given by ( 18 ) . If the origi- 
nal network of springs homogenizes to an effective elasticity tensor field C(x) 
then the new network of torque springs homogenizes to an effective elasticity 



tensor field C'(x) given by (11 1, assuming the transformation x'(x) only has 
variations on the macroscopic scale. In particular the stress field in the ho- 
mogenized network of torque springs is not be symmetric, and is influenced 
not just by the local strain, but also by the field of microrotations. 

There are some practical barriers to this homogenization. Suppose, for 
simplicity, that we are in two dimensions, that the original network consists 
of a triangular network of identical springs with bond length h under uniform 
loading so that the tension is the same in all springs, and that the transfor- 
mation is a rigid rotation x' = Rx where R^R = I. The displacement Uj of 
the nodes x^ is, up to a translation, that of uniform dilation, = axj. It 
follows that if i and j are adjacent nodes on the network, then — u^- scales 
in proportion to h. On the other hand, in order that the traction force per 
unit length on a line remains constant the tension T in each torque spring 
must also scale in proportion to h. Therefore the torque spring constant 
k' = Kmuj^ / (moj^ — 2K) must be essentially independent of h. Also we don't 
want the density of mass per unit area associated with the torque springs to 



10 



F. Guevara Vasquez, G.W. Milton, D. Onofrei and P. Seppecher 



be too large (otherwise gravitational forces would be very significant). This 
would be ensured if m scales as where /? > 2. Since 

fcW^ (21) 
2k' + TOw^ 

we see that K should also scale as h'^, and that 2K would be close moj^ 
when h is small. Thus each torque spring is very close to resonance. If this is 
satisfied at one frequency, it will not be satisfied at nearby frequencies. Thus 
the metamaterial is operational only within an extremely narrow band of 
frequencies. The situation is similar in three dimensions in a network having 
bond lengths of the order of h. Then — u^-, T, k' and m need to scale as 
, h, and , respectively, with /3 > 3 to avoid an infinite mass density 
in the limit /i — 0. (T must scale as h'^ to maintain a constant traction per 



unit area on a surface). Again K given by (21 ) must be close to mw^/2 when 
h is small. 

In three dimensions an alternative is to avoid the use of masses within each 
torque spring altogether. This can be achieved by pinning the internal nodes 
of the torque springs, where there would be masses (such as at the nodes ti 
and t2 in Fig.[T]), to a rigid lattice (designed in a way which avoids intersection 
with the springs inside the torque springs). Such a pinning corresponds to 
setting m — oo and each torque spring has then a spring constant k' — K 
which is independent of frequency. The resulting metamaterial is operational 
at all frequencies. Note that within the framework of linear elasticity each 
torque spring exerts a torque but not a net force on the underlying rigid 
lattice. If the rigid lattice (which might have only finite extent) itself is not 
pinned we require that the external forces on the metamaterial to be such 
that there is no net overall torque on the rigid lattice. 

A more serious concern is the validity of linear elasticity, at least using 
the torque spring designs proposed here. A characteristic feature of the de- 
signs involving masses is that the internal masses m do not move when the 
springs are translated, to first order in the displacement. This accounts for 
the balance of forces ^ = — y However the masses do move significantly 
if the terminals are translated a distance which is comparable to the size of 
the torque spring. Alternatively, if we pin the internal nodes of the torque 
springs, where there would be masses, to a rigid lattice then this restricts 
the motion of the torque spring terminals relative to the lattice. Clearly for 
the operation of the metamaterial the displacements must be small com- 
pared to h, assuming the size of each torque spring is of order h. When h 
is very small this severely limits the amplitude of waves propagating in the 
metamaterial for which linear elasticity applies. Thus the only metamaterials 
of the type described here that might possibly be of practical interest are 
those for which h is not too small. This is in contrast to honiogenization of a 
normal elastodynamic network where linear elasticity may apply when only 
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the displacement differences u- — u^-, between adjacent nodes i and j, are 
small compared to h. 



2 Active exterior cloaking in the quasistatic regime 

We show that for the Laplace equation, it is possible for a device to generate 
fields that cancel out the incident field in a region while not interfering with 
the incident field far away from the device. Our results generalize to the 
quasistatic (low frequency) regime. Thus any (non-resonant) object located 
inside the region where the fields are negligible interacts little with the fields 
and is for all practical purposes invisible. Here we relate the problem of 
designing a cloaking device to the classic problem of approximating a function 
with polynomials. Then we propose a cloak design that is based on a family 
of polynomials. We also show that our solution can be easily modified to give 
cloak objects while giving the illusion of another object (illusion optics as in 



2.1 Active exterior cloak design 

Following the ideas presented in [12] , we first state the requirements that the 
field generated by a device (source) needs to satisfy in order to cloak objects 
inside a predetermined region. Here we denote by i?(x, r) C the open ball 
of radius r > centered at x e . 

Let B{c,a) with a > and c e be the region where we want to 
hide objects (the cloaked region). The cloaking device is an active source 
(antenna) located (for simplicity) inside B{0,5) with S 1. Assuming a priori 
knowledge of the incident (probing) potential uqj we say that the device is an 
active exterior cloak for the region B{c,a) if the device generates a potential 
u such that 

i. The total potential u + Uq is very small in the cloaked region B{c, a). 

ii. The device potential u is very small outside B{0, R), for some large R> 0. 

Therefore, if the incoming (probing) field is known in advance, an active 
exterior cloak hides both itself and any (non-resonant) object placed in the 
region B{c, a). Indeed any object inside B{c, a) only interacts with very small 
fields and the device field is very small far away from the device. 

After a suitable rotation of axes, we may assume, without loss of generality, 
that c = (p, 0) with p > 0. As in [T^, we require the following conditions in 
our cloak design, 
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p > a + S, the active device is outside the reEion B(c,a), and 

e, \ , J, (22) 

R > a + p, the cloaking effect is observed in the far field. 



2.2 The conductivity equation 

Next, in the spirit of [TH [23] we give a more rigorous formulation of the 
exterior cloaking problem for the two-dimensional conductivity equation and 
prove its feasibility. The results extend easily to the quasistatic regime. 



r 

/ -0(0,5) 


B(c, a) \ 


c 

L 





w = l/z 



z = 1/w 



/®(o,i) 

— P-^ 









Fig. 3 The regions of Theorem[l](left) and their transforms under the Kelvin or inversion 
transformation (right). The device field is harmonic everywhere except in the gray areas, 
is close to minus the incident field in the red region and is close to zero in the green region. 



Theorem 1. Let a, c, R and S satisfy (22), then for any e > and any 

p2 



harmonic potential uq, there exists a function go 
u : — >■ K, satisfying 



and a potential 



Au^O, mM^\B(0,^), 

u — go, on dB{Q, 5), 
\u\ <e in R^\B{0,R), 



(23) 



- Uq\ < e in B{c, a) 



Proof. By applying the inversion (or Kelvin) transformation w = l/z, the 
geometry of problem ( 23 ) transforms as follows, 

• R'^\ B{0, S) transforms to B{0,l/d), 

• M2 \ 5(0, R) transforms to B{0, 1/R), 

• B{c,a) transforms to B{c*,a), with 
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a 



— I ' 



(^,0), and P 



The different regions and their transforms are ihustrated in Fig. [3j Thus the 
problem ( 23 1 is equivalent to finding 50 smd u such that 



AS = 0, in 5(0,1/(5), 
" = 5o, on 95(0, l/J), 



|m| < e, in B(0,l/i?), 



(24) 



^ \u + uqI < e, in i3(c*, a). 



Relating to the functions go and uq from (23), we get 5o(^) = 50(1/2) and 
uq{z) = uo(l/-2^)j so that uq is harmonic in the whole space except the origin. 
Next, we observe that the inversion transforms the necessary conditions (22 1 
to 



< /3 — a, the two balls 5(0, 1/i?) and B{c* , a) do not touch, 
1 



(25) 



P + a < -, the two balls 5(0, 1/(5) and 5(c*,a) do not touch. 




Let Uq be the analytic extension of uq in 5(c*, a), obtained with the harmonic 
conjugate such that up is the real part of Uq. Because of analyticity of Uq, 
we can approximate Ug with a polynomial Qq (e.g. by truncating the series 
expansion of Uq) such that 



M < 2 in B{c*,a). 



(26) 



This immediately yields the approximation for uq 



Wo - gol < 2 in B{c*,a), 



(27) 



where qo = 5R((3o), i-c, the real part of Qq. Since Uq can be approximated 
arbitrarily well by a polynomial, it is enough to consider (24 1 when uq is the 
real part of a polynomial, i.e. 



Au^O, in 5(0, 1/(5), 
u = go, on 95(0,1/(5), 



\u\ < e, in 5(0, 1/i?), 



(28) 



[ \u + qo\ < e/2, in 5(c*,a)). 



In other words, problem (28 1 is equivalent to finding a function u, harmonic 
inside 5(0, 1/(5) that approximates qq well inside B{c*^a) but is practically 
zero in 5(0, 1/R). 
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Let us now recall a classic result in harmonic approximation theory due 
to Walsh (see [Tl|, page 8). 

Lemma 1 (Walsh). Let K be a compact set in such that M.'^ \ K is 
connected. Then for each function w, harmonic on an open set containing K , 
and for each d > 0, there is a harmonic polynomial q such that \w — q\ < d 
on K . 

Walsh's lemma implies the existence of a harmonic solution to problem 



(28). Indeed, from the design requirements (25) there exists < ^ <C 1 such 
that 

-^+e</3-a-e (29) 

Then applying Lemma [l] with K = B{0, l/R) U B{c* , a), we obtain that for 
an arbitrary small parameter < <C 1 and for the function w satisfying 

.=r (30) 

^ - inB(c*,a + e), 
there exists a harmonic polynomial q such that Ig— wj < d on K. We conclude 



that there exists a harmonic solution to problem ( 28 ) , which implies the 
statement of Theorem [T] □ 



2.3 Explicit polynomial solution in the 
zero frequency regime 

Although mathematically rigorous, the existence result of Theorem [l] (which 
follows from Walsh's lemma) does not give an explicit expression for the 
required potential at the active device (antenna). In [T^ (see also [23 ) we 



give a polynomial solution to problem ( 24 1 . Unfortunately the radius a of the 



cloaked region in the polynomial solution of [191 ES] is limited by the distance 
from the origin p according to a < (2 + 2y/2)~^p. Thus in (THl [53] we can 
only cloak large objects if they are sufficiently far from the origin. Here we 
state a conjecture that extends our previous results [T^ [53] and that gives 
more freedom on the choice of the cloaked region location and size. This is 
supported by numerical evidence (see Figs. |4] and [5| . 

Conjecture 1. Let c* — (/?, 0) be as in the proof of Theorem[l] For any L > 0, 
any disk Si in the connected component containing the origin of the set 

Dp,L = |z e C, |z ~ < (1^1^^^} , (31) 



any disk >5'2 in the connected component of the set Dj^x containing the point 
c* and any e > 0, there exists two positive integers s and n such that \s/n — 
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L| < e and the polynomial Pn,. 



defined by 



P„,,(z) = 1 - 



s n — 1 

E 



s + j 
j 



satisfies 



11 < e on dSi and |P„ s| < e on 95*2. 



(32) 



(33) 



Moreover the approximation property ( 33 ) is not satisfied when either Si or 
1S2 is not contained in D^ l. 



Remark 4- To see why we expect that the polynomial P„ ^ satisfies (33 1, no- 
tice that ^ = /3 is a root of multiplicity s of the polynomial Pn^s- From the 
Taylor expansion of P^ g around z = f3, we can expect that P„_s « in a 
sufficiently small disk around /3. (The symbol ~ denotes an approximation 
with respect to the supremum norm.) Now the function 



E 



s + j 
j 



1 



has a root of multiplicity n at z = 0, i.e. ^''^•'(O) = for fc = 0, . . . , n — 1. This 
is because the sum in the definition of g{z) corresponds to the first n terms 
in the Taylor expansion around z = of (1 — z/ j3)^^ . Thus by Leibniz rule, 
z = is a root of multiplicity n of the polynomial P„ ^ — 1 = (1 — z//3yg{z), 
and we can expect Pn.s ~ 1 in a sufficiently small disk around the origin. This 
suggests an alternative definition of Pn.s as the unique Hermite interpolation 
polynomial (see e.g. |lHj) satisfying: 

-Pn,.(0) = 1, 

P«(0) = 0, forfc = l,...,n-l, 
pW(/3) = 0, forfc = 0,...,s-l. 

Remark 5. To motivate our belief that the region is the region of con- 

vergence of Pn,s as n — >■ 00 and s — >■ 00 with s/n ^ L, consider the special 
case where L > is an integer and s = nL. Then the last term in the sum 



(321 defining Pn,nL{z) is 



n(L + 1) - 
n — 1 



which diverges outside of Dpx as n — > 00 because 
"n(L +1)-1\^ (L + l^^+i 



Li 



, as n cx), 
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which follows from Stirling's formula, see e.g. f4T, §5.11]. Therefore the region 
of divergence of Pn,Ln contains the complement of Dpj^. 



For some polynomial Qo we can deduce from ( 33 ) that 



(34) 



QaPn,s - Qo ~ on dSi and QaPn,s - Qa ~ -Qa on 882, 



Thus, by construction, the real part of = QaPn^s — Qo is a solution of ([28j) 
with 5*1 = 5(0, 1/i?) and 6*2 ~ B{c* ,a). In the particular case when n = s 
(i.e. L ~ 1), Conjecture [1] was proved in [23], see also [19]. 




(a) n 



o 



(b) n = 5, s = 25 



Fig. 4 Contour plot for the polynomial Pn.s with /3 = 1. The solid white line is the level- 
set |P„ 3(2)1 = 10~^, thus the cloaked region could be any disk inside this level-set. The 
dashed white line is the level-set \Pn^s{z) — 1| = 10"'^. The device field is small in any 
circle inside this level-set. The red curve is the boundary of Dp i^, the conjectured region 
of convergence of Pn,B as both n — > 00 and s — > 00 with s/n = 1 and 5, respectively. (This 
is proved in the case n = s in 1231 .) The color scale is logarithmic from 0.01 (dark blue) to 
100 (dark red), with light green representing 1. 



In Fig. |4]we present a contour plot of the polynomial P„_s when /? = 1 for 
different values of n and s. The region bounded by the peanut shaped red 
curve represents the conjectured domain of convergence of the functions 
Pn,s when n — )■ 00, s — )■ 00 and s/n — >■ L. In the left side of D^j,, Pn,s 
is conjectured to converge to one, while in the right side of -D^,l, Pn,s is 
conjectured to converge to zero. The area within the solid white circle on 
the right represents the region to be cloaked and the area within the dashed 
white circle in the left represents the location of the observer. We now present 
an active cloak design based on Conjecture [l] 

Remark 6. Let uq be an a priori determined incoming harmonic potential. 
Let n and s be such that the real part of = QoPn,s — Qa is a. solution 
of (28) (recall Qq is a polynomial approximation of Uq in B{c*,a)). Let S 
be a bounded region in the complex plane compactly including the two disks 
Si = -6(0, j^) and 52 = B{c*,a). Then, the cloaking strategy we propose 
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10| 




(a) n = s = 15 




(b) n = 5, s = 25 



Fig. 5 Real part of the total field with the cloaking device active, incident field uo{x, y) = x 
and /3 = 1. The solid white, dashed white and red lines are the inversion (Kelvin) transforms 
of their counterparts in Fig. |4] The small black scatterer inside the cloaked region is an 
almost resonant disk centered at (p, 0) with radius r and dielectric constant e given by: (a) 
p = 1.1, r = 0.2 and e = -0.99; (b) p = 1.7, r = 0.9 and t = -0.998. In (a) the color scale 
is linear from -10 (dark blue) to 10 (dark red). In (b) the scale is linear from -35 to 35. 



consists of an active device (antenna) located inside B{0,d) and capable of 
generating a potential equal to the real part of W{^) on the set G C, -j G 
dS}. By (27), the total potential in the original physical configuration (the 
field from the antenna plus uq) is well approximated by the real part of 
{W + Qo){l/z) which ensures an almost zero field region in B{c,a) with 
negligible perturbations on the field outside B{0, R). 

Fig. [5] illustrates how the cloaking device (represented by the solid black 
curve) works after applying the back-inversion to the configurations presented 
in Fig. [4] Here the incident field is uo{x,y) — x and the objects we want to 
hide are almost resonant disks. Clearly, the active device generates the nec- 
essary field to cancel the field in the cloaked region while having a very small 
effect in the far field (outside the white dashed circle). With a polynomial 
of the same degree, when s = 5n (Fig. [5jb)) we can hide an object roughly 
four times larger than when s = n (Fig. [5ja)). Thus using the polynomials 
Pn,s allows us to cloak large objects without restrictions on the distance from 
the origin /? as was the case in [121 123] • The disadvantage is that cloaking 
is enforced on dB{0, R) (dotted white line in Fig. [5| with a larger R in the 
asymmetric L > 1 case than in the symmetric L = 1 case. For example to get 
a device field such that |m| < fO~^, R needs to be roughly five times larger 
when s — 5n (Fig. [sjjb)) than when s = n (Fig. [sja)). 
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2.4 Extensions and applications 

We now extend the previous results to the case of an incoming field having 
sources in \ B{0,R). 

Remark 7. The case studied in Theorem [T] (with an explicit solution in Con- 
jecture [l]) corresponds to an incoming field uq generated by a source located 
at infinity. The more general case corresponding to an incoming field having 
sources in \ i3(0, R) can be treated similarly. Indeed, the problem remains 



to find go and u satisfying (23), or equivalently go and u satisfying (24) where 
inside -6(0, l/R), uo{z) ~ uq{1/z) is harmonic. We can still approximate its 
analytic extension Uq by a polynomial in B(c*,a) and the proof goes as in 
Theorem [T] 

Although our main focus here is cloaking, the same ideas can be applied 
to illusion optics, where one wants to conceal an object by imitating the 
response (scattering) of a completely different object. 

Remark 8. Let ui be the response of an object we wish to imitate, i.e. an 
arbitrary potential harmonic in a set Di C such that \ 5(0, R) Di. 
Assuming the same notations as before, for any (known a priori) probing field 
Uq, harmonic in M^, there exists a function g E C{dB{0, S)) so that the field 
u generated by the active device (antenna) located in B{Q,6) satisfies: 

i. The total field u + uq is very small in the cloaked region B{c, a). 

ii. The device field u is close to Ui in \ B{0, R). 

Remark [s] follows from the inversion (Kelvin) transform and Lemma [I] by 
using an argument similar to the proof of Theorem [T] Using ideas similar to 
those in Remark [7j the result of Rem ark [8] can be generalized to the case of 
an incoming field with sources in \ B{0, R). 

To illustrate Remark [8] assume that the field ui is chosen to be the response 
field of an inhomogeneity I when probed with the incident field uq. Then 
Remark |8] means that when probing with the field uq, an observer located 
in the far field detects the inhomogeneity I regardless of the inclusion inside 
B{c,a) and without detecting the active illusion device. This creates the 
illusion that the object inside B(c, a) is the inhomogeneity I. 



3 Active exterior cloaking for the Helmholtz equation in 
three dimensions 

Previously in [1^1 HOj we designed cloaking devices generating fields close 
to minus the incident field in the region to be cloaked and vanishing far 
away from the devices. Miller [32] proposed an active cloak based on Green's 
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identities: a single and double layer potential is applied to the boundary of 
the cloaked region to cancel out the incident field inside the cloaked region, 
while not radiating waves. The idea of using Green's identities to cancel out 
waves in a region is well known in acoustics (see e.g. [131 1311 lU). Jessel 
and Mangiante [23] showed that it is possible to achieve a similar effect to 
Green's identities (and thus cloaking) by replacing the single and double layer 
potentials on a surface by a source distribution in a neighborhood of the 
surface. What makes our approach different is that the cloaking devices are 
multipolar sources exterior to the cloaked region and thus do not completely 
enclose the cloaked region. In [THl [2D] the cloaking devices are determined 
by solving numerically a least-squares problem with linear constraints. Our 
cloaking approach easily generalizes to several frequencies [20j but requires a 
priori knowledge of the incident field. Zheng, Xiao, Lai and Chan [ST] used the 
same principle to achieve illusion optics |27| with active devices, i.e. making 
an object appear as another one. Then in [55| we showed Green's identity 
can be used to design devices which can cloak or give the illusion of another 
object, i.e. achieving an effect similar to the active devices in [19l [20l IST] . 
The single and double layer potential needed to reproduce a smooth field 
inside a region while being zero outside is given by Green's identity and can 
be replaced by a few multipolar sources using addition formulas for spherical 
outgoing waves. If in addition we want to imitate the scattered field from an 
object as in [51], a similar procedure applies. 

The active cloaking devices we designed in [l9l |2Ql |22] are two dimen- 
sional. Here we extend the result in [52] to the Helmholtz equation in three 
dimensions. The wave pressure field u(x) solves the Helmholtz equation, 

Ait-l-fc^M^O, forxeM^ 

where k — 27r/A is the wavenumber, A = i-nc/ui is the wavelength, c is 
the wave propagation speed (assumed to be constant) and lo is the angular 
frequency. Recall for future reference that the radiating Green's function for 
the Helmholtz equation in three dimensions is 

exp[zfc|x-y|] 

Gx,y^— — — 35 

47r|x~y| 

Another underlying assumption is that the frequency w is not a resonant 
frequency of the scatterer we wish to hide. 



3.1 Green's formula cloak 

As pointed out by Miller [32] it is possible to cloak an object inside a bounded 
region D G M'' from an incident wave (probing field) Ui by generating a 
cloaking device field using monopole and dipolc sources (single and double 



20 



F. Guevara Vasquez, G.W. Milton, D. Onofrei and P. Seppecher 



layer potential) on dD. The device field Ud can be defined using Green's 
formula 

Md(x) = / ASy {-(n(y) • VyU,(y))G(x, y) + u,(y)n(y) • VyG(x, y)} 

JdD 

otherwise, 

(36) 

so that the total field Ui + is a solution to Helmholtz equation for x ^ dD 
that vanishes inside D while being indistinguishable form Ui outside D. Since 
the waves reaching a scatterer inside the cloaked region D are practically zero, 
the resulting scattered field is also practically zero. For clarity we assume the 
region Z? is a polyhedron. The arguments we give here can be easily modified 




for other domains with Lipschitz boundary, as Green's identity (36) is valid 
for these domains [11]. 



Remark 9. The Green representation formula ( 36 ) requires that Ui be a C 
solution to the Helmholtz equation inside D. A similar identity holds when 
Ui is a radiating solution to the Helmholtz equation outside D. In this 
case, the device field Ud vanishes inside D and is identical to —Ui outside 
D. The exterior cloak we present here can in principle be used to conceal a 
known active source and possibly accompanying scatterers inside D. If the 
radiating wave Ui is taken to be the scattered field from a known object, the 
same principle can be used for illusion optics [UI [5T] . 



3.2 Active exterior cloak 



The main idea here is to achieve a similar effect to the Green's identity cloak 
but without completely surrounding the cloaked region by monopoles and 
dipoles on dD. We "open the cloak" by replacing the single and double layer 
potential on each face dDi of dD by a corresponding multipolar device located 
at some point x/. Each device produces a linear combination of outgoing 
spherical waves of the form 

«<iW = EE E 6;.«.™C'(x-xO, (37) 

l—l n— m— — n 

where Udev is the number of devices (or faces oidD) and V,"\x) is a radiating, 
spherical wave defined for x 7^ by 

yr(x) = /iW(fc|x|)r„"(x). 
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Here h^n \t) is a spherical Hankel function of the first kind (see e.g. [13 
§10.47]) and Y™{x.) is a spherical harmonic evaluated at the point x = x/ |x| 
of the unit sphere 5(0, 1). In spherical coordinates, the spherical harmonics 
we use are defined as in [51 §2.3] by 

>T(M)^^?^|^el(cos.)e'-, (38, 

where the elevation angle is 9 G [0,7r] and the azimuth angle is </) <E [0, 27r]. 
Here /^"'(t) are the associated Legendre functions 

defined for n — 0, 1, 2, . . . and m — 0, 1, . . . , rt in terms of the Legendre 
polynomials P„ of degree n with normalization ^"^(1) — 1- The definition 
(38) ensures that the spherical harmonics have unit ^^(^(0, 1)) norm. 

The main tool to replace the fields generated by a face is the addition 
formula (see e.g. Theorem 2.10 in [5]) 

CO n 

G(x,y)=zA:^ ^ y„™(x)c7™(y) (39) 

n— m— — n 

which means we can mimic a point source located at y by a multipolar source 
located at the origin. The coefficients in the multipolar expansion are values 
of entire spherical waves 

c/r(x)=j„(fc|x|)rr(X), 

where jn(i) are spherical Bcssel functions [HJ §10.47]. The series in the mul- 
tipolar expansion (39) converges uniformly on compact sets of |x| > |y|. 
We are now ready to state the main result of this section. 



Theorem 2. Multipolar sources located at the points x; ^ dD, Z = 1, . . . , n^, 

can be used to reproduce the Green's formula cloak outside of the region 

A= [J B (x,, sup |y-Xi| 1 , 

where B{:x.,r) is the closed ball of radius r centered at x. The coefficients 
bi^n,m (37 1 such that u[l^*\x) = Ud(x) for x ^ A are 

bi,n,m ^ik ( dS-yf (-n(y) • VyU,(y))C/;»(y - x,) 

JdD, (40) 



+ 7.,(y)n(y).VyC/™(y-x,)}. 
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Moreover the convergence of (37) is uniform on compact sets outside A. 



Proof. Splitting the integral in ( 36 ) into integrals over each of the faces dDi 
of the polyhedron dD and applying the addition theorem ( 39 1 with center at 
the corresponding X; we obtain: 



n—0 m——n 
oo n 



n—0 m— — n 



(41) 

The result (40) follows for x ^ A by switching the order of the sum and 
the integral in (41). For the first term in the integrand of (|4T|), this switch 



is justified by the uniform convergence of the series (39) (for all devices) in 
compact sets outside of A. 

For the second term in the integrand of (41 ), we shall show that the series 
converges uniformly on compact sets outside A, so it is also valid to switch 
the integral and the series in (41 ). To see the uniform convergence, it is useful 
to split the products y^'(x — Xi)VyC/™(y — x;) into two terms corresponding 
to the two terms in the gradient 

Vyf/„"(y) = fcyj;(fc|y|)^r(y) +j«(fc|y|)^^^^J^3^(vrr)(y) ,,,, 

|y| (42) 

= &{y) +5S„(y), 

where / is the 3x3 identity matrix. 

For the series involving the first term in the gradient ( 42 ) we bound with 
the triangle and Cauchy-Schwarz inequalities: 



m— — n 



^ K"(x-x05il(y-xO 

n 

/.«(A:|x-x,|)j;(fc|y-xn)| ( ^ 



< k 



Using the summation theorem for spherical harmonics (see e.g. Theorem 2.8 
in[5]) 



m— — n 



2n + l 
in 



for any y € ^(0, 1) and n = 0, 1, 



(43) 



we get the estimate: 
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m— — n 



<ft kW(fc|x-x,|)j;(fc|y-xz|) 



2n+ 1 
47r 



(44) 



|y-x/| 



in+l I ' 



for large n — >■ oo. The last equality comes from the asymptotic expansion of 
Bessel functions for fixed t > and large order n, (see e.g. [41', §10.19]) 



and 



/.W(t) 1). 



For the series involving the second term in the gradient ( 42 1 we bound the 
sums 



^ yr(x-x05il(y-xO 



m— — n 



< 2 



ftW(fc|x-x,|) 



j„(fc|y-X;|) 



|y-Xi| 



Using the summation theorem for spherical harmonics ( 43 ) and their gradi- 
ents (see e.g. (6.56) in ^8j), 



m— — n 

we get the asymptotic 



n{n + l){2n+ 1) 



An 



for any y e S'(0,1), (45) 



E K"(x-xOffil(y-Xi; 



m— — n 



< 2 



/iW(fclx-xn) 



i«(fc|y-x;|) 



2n + ly f n{n+ l)(2n + 1) 



47r 



47r 



(46) 



|y-x/i 



x-Xi 



iri+l 



Here we have used that for t > fixed and as n oo, (see e.g. [HI §10.19] 

|j„(t)|=0(r) and hll\t) =0(i-"-i). 
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The estimates (44) and (461 give uniformly convergent majorants for the 
series in the second term of (41), since when x ^ A we have |y — x/| < 
for / = !,..., rirfp,,. The proof is now completed. □ 



3.3 A family of exterior cloaks with four devices 




Fig. 6 The configuration for the tetrahedron based cloak of Sect. 3.3 The distance in 
red is the radius a of the circumsphere to the tetrahedron D. The distance in black is the 
distance 5 from the origin to a device. The distance r{a, S) (in green) is the distance from 
a device to the closest vertex of D. The exterior surface of the region A of Theorem [2] IS 
in grey and has been cut to reveal the cloaked region D \ A in red. The four devices are 
shown with stars. 



Nothing in Theorem [2] guarantees that the cloaked region D \ A is non- 
empty. We show here how to construct a family of cloaks with non-empty 
D\A based on Green's identities applied to a regular tetrahedron D. We also 
determine what is the position of the devices that gives the largest cloaked 
region within this family. 

Consider a regular tetrahedron with circumsphere S'(0, a) and vertices 
ai, . . . ,a4. We locate the devices Xi, . . . ,X4 on 5(0,(5), with S > a, such that 
X; replaces the face opposite to vertex a;, that is X; and a; are on opposite 
sides of the plane formed by the face of the tetrahedron not containing a;. 
For simplicity we also require that x; — a; is normal to this plane. The con- 
figuration is sketched in Fig. [6j Simple geometric arguments show that the 
radii of the balls that define the region A are all equal to 
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r{'7,d)= \ [a--] +-6'] . (47) 

Moreover the radius of the largest sphere fitting inside the cloaked region is 

r,sia,6)=S-r{a,S). (48) 

For fixed S, the largest possible cloaked region is obtained when a = 6/3 
which corresponds to the case when every triplet of balls in the definition of 
region A touch at a vertex a; of the tetrahedron. Thus for fixed 5, the largest 
sphere we can fit inside the cloaked region has radius, 



0.057(5. (49) 




3.4 Numerical experiments 

We report in Fi g. [7 simulations of this cloaking method with the setup de- 



scribed in Sect. 3.3 



The incident field we take is the plane wave iti(x) 
exp[iA;k • x] with direction vector k = [1,1,1]/a/3. We first compute the 
device field of Theorem [2] by truncating the sum in n of (37 1 to n < N. 



Throughout our numerical experiments we determine TV with the heuristic 
(found by numerical experimentation) 

7V((5) = [f.SH], (50) 

where \x~\ is the smallest integer larger than or equal to x. The integrals in 



( 40 ) were evaluated with a simple quadrature rule that is exact for piecewise 
linear functions on a uniform triangulation of the faces of the tetrahedron 
D, we chose the number of quadrature points so that there are at least eight 
points per wavelength. The scattered field by a ball was computed by first 
evaluating the incident field (or device field depending on the case) on a grid 
with equal number of points in (j) and 9 and then finding its first few spherical 
harmonic decomposition coefficients using the sampling theorem [TO] . 

As can be seen in the first row of Fig. [7] the device field Ud is virtually zero 
far from A while being close to the incident field in the cloaked region D\A. 
In the second and third rows of Fig. [T] we display the total field in the presence 
of a sound-soft (homogeneous Dirichlet boundary condition) ball centered at 
the origin and of radius 3r*g{d) (i.e. a larger scatterer than what we expected 



from Sect. 3.3). The scattered field from the ball reveals the ball's position 
when the devices are inactive (third row). The scattered field is essentially 
suppressed when the cloaking devices are active (second row), as the field is 
indistinguishable from a plane wave far from A. 
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Since as i ^ 0, h^n\t) = ©(i""-!) (see e.g. PH §10.52]), we expect the 
device field to blow up as we get close to the device locations X;. This 
blow up corresponds to the "urchins" in the first and second rows of Fig. [7| 



where even with the truncation of the series ( 37 1 , we observe very large wave 
amplitudes which would be hard to realize in practice. Fortunately we can 
enclose the regions with very large fields by a surface and apply Green's 



formula (36 1 to replace these large fields by (hopefully) more manageable 
single and double layer potentials on the surface of some "extended" cloaking 
devices. 

We illustrate these "extended" devices in Fig. [8] where we display the level 
sets where the device field amplitude is 5 (or 100) times the amplitude of the 
incident field. At least for the particular configuration {S — 6A) considered 
in Fig. |8] these surfaces resemble spheres surrounding each device location 
X/. The "extended" devices still leave the cloaked region (in red in Fig. [8| 
communicating (connected) with the background medium. This is why we 
call our cloaking method "exterior cloaking" . 

We also consider the extended devices for larger values of S in Fig.|9] Here 
we look at the cross-section of the extended devices on 5(0, cr), which in the 
construction of Sect. |3.3| is the circumsphere to the tetrahedron D. In the 
optimal case S — 3a, the predicted cloaked region D \ A and the exterior 
K'^ \ A meet on 5'(0, a) at the vertices of the tetrahedron D. We see that the 
extended devices (in black in Fig. [o]) grow as S increases, and leave gorges 
communicating the cloaked region with the exterior. The centers of the gorges 
appear to agree with the vertices of the tetrahedron D. The percentage area 
of S{0, a) that is not covered by the cross-section of the extended devices on 



5(0, cr) is also quantified in Fig. 10 h). Since the relative area of the openings 



appears to decrease monotonically with S/X, Fig. 10 b) suggests the gorges 
close for large enough 5/X. Further investigation is needed to find out whether 
the shrinking openings in the cloak is due to our choice of N with heuristic 



(50). 



Finally we give in Fig. 10 ^a) quantitative measures of the cloak perfor- 
mance for different values of S. These measures show that the device field is 
close to minus the incident field inside the cloaked region and that it is very 
small outside of the cloaked region. 
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z = —2a z = —a 2 = z = a z = 2a 




Fig. 7 Constant z slices of the real part of different fields, for the optimal case 5 = 3a 
and with <5 = 6A. The first row shows the device field which is close to zero far from 
the devices and close to —Ui in a small region close to the origin. The second and third 
rows show the total field when the devices are active and inactive in the presence of a 
scatterer. The scatterer is a sound-soft ball centered at the origin and of radius 3r*g{5). 
Even though this ball is not completely contained inside the tetrahedron D, the scattered 
field is greatly suppressed when the devices are active, making the ball harder to detect 
far from the devices. The color scale is linear from -1 (dark blue) to 1 (dark red) and each 
box is lOA X lOA, with the z— axis at the center. 




(a) \ua\ = 100 (b) \ua\ = 5 



Fig. 8 Contours of \u^\ (gray) and \uj^ + Ui\ = 10^^ (red). Here the vector (0,0,1) is 
perpendicular to plane of the page. By Green's identity it is possible to replace the large 
fields inside the gray surfaces by a single and double layer potential at the gray surfaces. 
These "extended devices" need only to generate fields that are at most the fields on the 
contours that we plot and they cloak the red region without completely surrounding it. 
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(5 = 6A 5 = 12A 




5 = ISA S = 24A 



Fig. 9 Cross-section of level set \u^\ > 10^ (black) and of the region A (shades of gray) 
on the sphere |x| = cr for the optimal rr = 5/3. Here we used the equal area MoUweide 
projection (see e.g. |12|). In the optimal case, each triplet out of the four balls forming A 
meets at a single point which is a vertex of the tetrahedron D. Note that for the cases in 
the first row there are four distinct extended devices. The leftmost and rightmost spots 
correspond to one single device split in two by the projection. 




6 12 18 24 6 12 18 24 



(a: <5/A) (b: <5/A) 

Fig. 10 (a) Cloak performance. In red: \\ui + where the norm is the 

L-^(5(0, r*g((5)) norm, which measures how well we approximate the incident field inside 
the cloaked region. In blue: where the norm is the L^(S(0,2<5)) norm, which 

measures how small is the device field far away from the devices, (b) Percentage of the area 
outside the cross-section of the extended devices on the sphere S{0,a = 5/3) for different 
values of (5. 
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